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Abstract 

Tripartite three mode Gaussian state undergoes parametric amplification and amplitude damping as well 
as thermal noise is studied. In the case of a state totally symmetrically interacting with the environment, 
the time dependent correlation matrix of the state in evolution is given. The conditions for fully separability 
and fully entanglement of the final tripartite three mode Gaussian state are worked out. 
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1 Introduction 

Quantum entanglement of continuous variables (CV)[T] |^ has received much attention recently, as shown by the 
spectacular implementations of deterministic teleportation schemes [3] [1] [5] [1] [7j, quantum key distribution 
protocols [8], entanglement swapping f5] [9] [10], dense coding [TT], quantum state storage [121 and quantum 
computation [13j processes in quantum optical settings . Quantum information with CV in general is mainly 
concerned with the family of Gaussian states, since these comprise essentially all the experimentally realizable 
CV states. A practical advantage of CV systems is the relative ease with which entangled states can be 
generated in the laboratory \jt\ . Former works are mainly on the entanglement of bipartite system. The study 
of CV multipartite entanglement which was initiated in [B] [S], where a scheme was suggested to create pure 
CV N-party entanglement using squeezed light and N-1 beam splitters. In the practical situation, such a pure 
multipartite entanglement state will evolve to a mixed state, due to the decoherence. In all the multipartite CV 
states, tri-mode entangled state is the simplest one, and a complete classification of tri-mode entanglement was 
obtained, directly computable criterion that allows to determine to which class a given state belongs[n] [T5] . 
We in this paper will investigate the separability of tripartite Gaussian state in presence of amplitude damping, 
parametric amplification and noise which are symmetric among the modes, based on our former works on the 
corresponding problem of bipartite CV systems [TB] [T7] . 

2 Time evolution of characteristic function 

The density matrix obeys the following master equation [TH] [IH] [10] 

+ (1) 

with the quadratic Hamiltonian 

H ^ h'^^-{r]jka]al - T^*kajak) (2) 

where 77 is a complex symmetric matrix. In the single-mode case, this Hamiltonian describes two-photon 
downconversion from an undepleted (classical) pump [20] ■ The full multi-mode model describes quasi-particle 
excitation in a BEC within the Bogoliubov approximation [21j . This item represents the parametric amplifier. 
While the amplitude damping is described by £, 
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where the Lindblad super-operator is defined as L\o\p = lopo^ — o^op — po^o. 

We now transform the density operator master equation to the difi^usion equation of the characteristic 
function. Any quantum state can be equivalently specified by its characteristic function. Every operator A G 
Biji) is completely determined by its characteristic function XA '■= tr\AD{p^ [52], where V{p) — exp{pa^ — p*a) 
is the displacement operator, with p, = [pi, p2, • • • , Ps] , a = [ai, 02, • • • , fls]"^ and the total number of modes is 
s. It follows that A may be written in terms of XA as [23]: A = J[Yli ^v^]x.a(m)2-'(— A*)- The density matrix p 
can be expressed with its characteristic function X- X = ^''[p^(m)] • Multiplying 'D(p) to the master equation 
then taking trace, the master equation of density operator will be transformed to the diffusion equation of the 
characteristic function. It should be noted that the complex parameters pj are not a function of time, thus 
^ = tr[^'D{p)], the parametric amplification part in the form of characteristic function will be [20j 



]k jk ^'^ 

The master equation can be transformed to the diffusion equation of the characteristic function, it is 

dx , 9x 9x 

-Ql = -UW3Q^^+V,kf^3Q^) (4) 

jk ^ " 

-^Er.{lM,l^ + (2^. + i)l/^.f)x}- 

Where we denote pj as \pj\e^^\ and we should take care about that the variables are pj and p* in the 
amplification while they are \pj\ ,6j in the damping. 

3 The parametric amplifier and the ampUtude damping 

Suppose the solution to the diffusion equation is 

x(m,M*,0 = x(^,^*,0)exp{i(^,-^*)(^ ^ 

%-f^n( « f* ){p*,-pf}, (5) 



where v = pM + p*N, with M and N being time varying matrices, a and /3 are constant matrices and = a, 
(5 = 0^. Then M and N are the solutions of the following matrix equations [H] 

where T = d^ag{Tl,T2^ ■ ■ ■ ,1"^}. While a and (3 are the solution of the following matrix equations 

2ria + 2a*r^-Tl3-l3T + Tw + wV = 0, (8) 
ra + ar-2?7*/?-2/3*r;-r(2n+l)-(2n+l)r = 0. (9) 

with n = diag{ni,n2, • ■ ■ ,ns}. The constant matrices a and (3 can be worked out as the solution of linear 
algebraic equations @ and (O. What left is to solve matrix equations (O and ([7|). There are two situations 
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that the equations are solvable. The first case is that all the modes undergo the same amplitude damping, that 
is Fi = r2 = • ■ ■ = Fs, thus F = Fil^. F commutes with any matrix. Equations ([6]) and ([7|) have solution 



M = e-^'^' cosh* {\r]\t), (10) 



_irtSinh(|?7lO. 



N = ILLL^rj, (U) 

\m 



where the matrix cosh and sinh functions are defined as [23] 

coshKI = / + + + , (12) 

The second case is that 77 is a real matrix while the amplitude damping can be different for each mode. The 
solution of equations ^ and ([7]) will be 

M - i[exp(-r7<-y)+exp(77t-y)], (13) 
N = i[exp(-7/<-y)~exp(7yt-y)]. (14) 

4 The separability criterion of tripartite Gaussian state 

The separability problem of the three mode gaussian state was perfectly solved [14J. The three mode gaussian 
states were classified as 5 different entangled classes |14j. But the states in this paper can be classified as 
3 different entangled classes: fully inseparable states, biseparable states, fully separable states. Following the 
notation of Ref. |14| . the CM acM now is replaced with 7, where 7 — 2acM, the partial transposition is denoted 
as Aj [j — 1,2,3 is one of the three parties, here each party consists one mode. If the canonical observables 
are arranged in the order of xi, a;2, X3,pi,p2,P3, one has Ai = diag{l, —1, 1, 1, 1, 1}, A2 = diag{l, 1, 1, —1, 1, 1} 
A3 — diag{l, 1,1,1, 1,-1}. The partially transposed CM will be 7^ — Aj'yAj. Denote 



i 

then criterion for fully inseparable state is 



-1 

1 



(15) 



7j ^ iJ, for aU j = 1,2,3. (16) 

Because of the symmetry of the 3ST state, the criterion can be simplified to for example ^ iJ ■ 

While for jj > iJ, ( j = 1,2,3), the state will be PPT tri-mode state, and it can be biseparable or fully 
separable. The criterion to distinguish the biseparable and fully separable states is p3] as follow. The CM 7 of 
PPT tri-mode state can be written of as 

where y4 is a 2 x 2 CM for the first mode, whereas S is a 4 x 4 CM for the other two modes. Define the matrices 
K and K as 

K = A-C—^—C'^, K = A-C — ^C^, (18) 
a - iJ B - iJ 

where J = J0 (— J) is the partially transposed J for two modes. 
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Then the condition of the PPT tri-mode state being fuUy separable is that if and only if there exists a point 
(y, z) e i?^ fulfilling the following inequality: 

min{trif , tr^} > 2x, (19) 
detK + 1 + L'^ { y, z )^ > x ■ tiK, (20) 
det^ + l + L^( y, z )^ > x-trK, (21) 

where x = + j/^ + z^, and L ~ {u~ w, 2Re(u)) , L = (u — w, 2Re(?J)) if K and K is written as 

Ineq. (fT9|) restricts {y, z) to a circular disk C, while Ineq. ((20|) and Ineq. (f2T|) describe eUipses £ and £' respectively. 
The existence of the point (y, z) then turn out to be the intersection of the ellipses £ and £' and the circular 
disk C. 

The intersection of the ellipses £ and £' is a range in the yz plane which is bounded by the elliptic curves 
d£ and d£' . In the cases considered in this paper, I{,e{v) = 0, 'R.e{v) = 0, the two elliptic curves d£ and d£' are 
described by 

AeiK +l + [u-w)y = {u + w)x, (23) 
det^ + 1 + (u - w)?; = {u + w)x, (24) 

i9£ and d£' are centered at y axis of yz plane. The intersection of d£ and d£' is the solution of these two 
equations as far as 

X > yTTy2. (25) 

Thus the condition of the existence of {d£) D {d£') is obtained. 

Two situations of the intersection of the ellipses £ and £' and the circular disk C should be considered. The 
first is {{d£) n {d£')) C C, in this case, fully separabihty condition is determined by ((^ and (^5]) . The 
second is {{d£) D {d£')) £ C, in this case, we should consider that if one of the tops of £0 £' is contained in C 
or not, the two tops are determined by equations (^5]) and ([^ separately by setting x = yjl + ip. 



5 The symmetric amplification and damping of tripartite Gaussian 
state 

We consider the totally symmetric amplification and amplitude damping among all three modes, that is 

1] ^ riol + 7]iS, (26) 

with the matrix S having its entries Sij = 1 for i j and Sij = for i = j {i, j = 1, 2, 3); Fi = r2 = Fa, which 
are redenoted as F, and ni — n2 — ns, which are redenoted as n. For simplicity, we only consider the case of real 
rj. The matrix rj can be written as rj = U C?7^^,with C — diag{irjo + 2771, 770 — ^71 , ^70 ^ '7i } = ^diag{C,o^ Ci; Ci}; a-nd 



U 



1 

71 



\/2 

V2 %/3 1 

^/2 -V3 1 



(27) 



Then 



cosh(77t) — [/(iia(ji{cosh(Coi')i cosh(^ii'), cosh((it')}C/ ^ 

= ^ [cosh(Coi') + 2 cosh(Cii')] + | [cosh(Cot') - cosh(Cit')] 



(28) 
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where t' = ^i, and there is the similar expression for sinh(77t). The Eqs. ([HI) and ^ can be simphfied in the 
symmetric case and the solution is 



a = {2n+l){ail + a2S), 
13 = (2n+l)(/3i/ + /32^), 



(29) 
(30) 



^ith ai = i(^ + ^), a, = i(^ - ^), /3i = i(^ + ^), = i(T 



Co 



From (O, the complex CM at any time t can be obtained in the real symmetric amplification and symmetric 
damping case, 



7c = 



M -N 
-N M 



[7c(0) 



a (3 
(3 a 



M 

-N M 



(31) 



where M — e ^i"* cosh(?7t), N — — e a'"* sinh(77i). The real parameter CM can be reduced from the complex 
CM and we have 



7 = n 
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(32) 



with n' = 2n+ 1, and a = iexp[2(Co - l)t'](l - y^) +f exp[2(Ci - !)<'](! - j^) +t(T4r + T^), 
6= iexp[-2(Co + !)<'](! -if^) +f exph2(Ci + l)i'] (1- ^f^) +2|^(^ + ^), c = i exp[2(Co - l)i'](l- j^) 
-i exp[2(Ci - mil - ^) +^(1^7 - T^): d = i exp[-2(Co + l)m ~ j^) -i exp[-2(Ci + l)i'](l - yfcr) 
+3"(iqr^ — Yqj^), where the vacuum initial state is assumed for simplicity. 

5.1 The fully separable conditions 

The condition ([^5]) of the existence of the intersection of the two ellipses can be simplified to 

- cd[{a - c)(6 + 2d) - 1][(6 - d){a + 2c) - 1] > 0. 
In the limit of t' — > oo, we obtain cd < for all the parameters Co f^nd Ci, denote 



a — a 



b' 
d' 



+ 2c = exp[2(Co - 1)<'](1 - + —— 

t — (,0 -L — <,0 

a - c = exp[2(Ci - 1)^'] (1 - 

ft 

fe + 2d = exp[-2(Co + l)t'](l - -) + 



1-Ci' 
n' 



6 - d = exp[-2(Ci + l)t'](l 



1 + Co 
n 



) + 



1 + Co' 
n' 



1 + Ci 1 + C; 



then the condition (l25ll turns out to be 



(a'd'- l)(c'5'- 1) > 0. 



(33) 

(34) 
(35) 
(36) 
(37) 

(38) 



When the amplification is weaker the the damping, that is, max{|Co| , |Ci|} < 1- In the limit of t' — s- oo, all 
the time dependent terms in a', 6', c' and d' are damped to zeros. The final CM tends to the residue CM, the 
condition of fully separability will be 



n'^ > (l + Co)(l-Ci), for Co > Ci, 
n" > (l-Co)(l + Ci), forCo<Ci, 



(39) 
(40) 
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Figure 1: The border of fully separable states and biseparable states, characterized by the noise, the ratio of 
single mode amplification to amplitude damping, the ratio of inter-mode amplification to amplitude damping. 
The noise is the average noise n. 

which can be rewritten as 

n'^ > (l-r;(, + r;;)(l + 77^ + 2r;;), for,?; >0, (41) 
> il-Vo-H)i^ + Va-v'i), ioTrj[<0, (42) 

respectively, with ?7o i = ^ ?yo,i- These are the fully separable conditions of the residue states and shown in 
Fig.l. Concerning the condition we consider the critical situation of tiq = (1 + Co)(l — Ci) for (q > ^i, the 
condition ifTH)) reduces simply to Co > Cii thus the state is fully separable in the critical situation. By physical 
consideration, a state with n' > tiq then is fully separable, because the state will be made more separable by 
adding the noise. The same conclusion is true for the case of Co < Ci- Hence the fully separability conditions 
are given by ([39]) and (|40|) in the situation of weak amplification. 

Moreover, we can prove that and ([3D]) are also fully separable conditions of the final state {t' ^ oo) for 
strong amplification situation. We provide the proof of one of the cases here, the other cases can be followed with 
the same method. The case we considered is Co > 1, |Ci| < I7 in the case we have a' +00, b' = Yf^,c' = 73^, 

d' = > 0, thus a'd' > 1 and the condition reduces to b'c' > 1, which leads to ([39l). The condition (fT9| 
can also be fulfilled for the critical state in the strong amplification situation. 

5.2 The biseparable conditions 

The biseparable condition ^A^iJ can be simplified to 

1 - {a'b' + 86'c' + 8a'd' + c'd')/% + a'b'c'd' > 0. (43) 

In the weak amplification case, max{|Co| , |Ci|} < 1- When t' 00, we have a' = jz^, b' — jj^, c' — jzr^i 
d' = . The solution of Inequality (|43|) is 

> 1 - ^ (Co + I6C0C1 + C?) + 3^ ICo - Ci I y^288 + Co'+34CoCi+C?, (44) 

which can also be written as 

n''>l- v'o' - vWi + Iv? + lWi\ ,J^2 + 4r,'^^ + Ar,',r,[^7r,[^ (45) 

In the strong amplification case, we first consider the situation of Co > li |Ci| < 1- When t' — s- 00, we have 
a' = exp[2(Co - l)t'](l + ^) + +00,6' = ^,c' = d' = Denote a' = n'ai, b' = n'b,, 

c' = n'ci , d' — n'di , then Inequality is 

1 — (ai6i + 85iCi + 8aidi + cidi)n' /9 + ai6iCidin'^ > 0. (46) 
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2io'r 2n,/r 

Figure 2: The border of biseparable states and fully entangled states, characterized by the noise, the ratio of 
single mode amplification to amplitude damping, the ratio of inter-mode amplification to amplitude damping. 
The noise is the average noise n. 

The left hand side contain the linear and cubic terms of n' implicitly due to ai. To solve the inequality, we just 
consider that ai does not depend on n' for the moment ( The solution to the equality part of (|46p is a kind of 
iteration solution), the solution to the quadratic inequality of n'^ can be obtained easily, at the limitation of 

oi — > +00, it is 

f2 1/1 8 , 

= (l-Ci)(l + Co)-^(l-Ci)(Co-Ci)- (47) 

The second situation we should consider is Co > 1; Ci < —1, we have ai — > +00, di +00, bi — j^:^, ci = 
at — > 00, thus the solution of Inequality reduces to 

n">-^ = l{l-Ci){l + Co)- (48) 
9oiCi 9 

Similar results can be obtained for other domains of the parameters Co and d , They are 

(49) 
(50) 
(51) 
(52) 

There are no restrictions to the n' in the domain of Ci > l,Ci > 1 and domain of Co < ^ l,Co < —1, which 
means the final states in these domains are always biseparable, in fact they are fully separable according to 
([39]) and (|40|). The synthesis of ((45|) and (l47 |) - (|52|) is the biseparable condition for all the parameters. We 
can transform the parameters to the single mode amplification parameter 770 and inter-mode amplification 
parameter 771 (appeared in the form of ryp = 2770/r, rj^ = 2r]i/T), and the biseparable condition is shown in Fig. 2. 
A comparison of the figures as well as the formula show that the fully separable and biseparable conditions are 
quite close with each other. The detail of difference of the two is displayed in Fig. 3 for the noiseless situation. 
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Figure 3: The noiseless situation, the soHd hne for the border of fully separable and biseparable states, the 
dashed line for the border of biseparable and fully entangled states. The difference of the two curve is very 
small, the detail of the difference is amplified by a factor of 100 and shown by the dashdot line. 

6 Conclusion 

We consider the effect of amplification, amplitude damping and thermal noise on tripartite three mode Gaussian 
state. The three modes are identical in the initial preparation and in the later interaction with the environment 
through amplification, damping and thermal noise. In such an assumption, the analytical expression of fully 
separable and biseparable conditions for the final tripartite three mode Gaussian state are obtained for both the 
weak and strong amplification situation. In the weak amplification, no further restriction on the initial state 
is required, besides the identical of the three modes. While in the strong amplification case, vacuum initial 
state is assumed for the simplicity of the description. The separability conditions are characterized by the 
ratio of the single mode amplification parameter with respect to the damping coefficient, the ratio of the inter- 
mode amplification parameter with respect to the damping coefficient, and the thermal noise. The biseparable 
condition and the fully separable condition arc very close with each other, make the domain of biseparable but 
not fully separable states quite small. While both the domains of the genuine entangled tripartite states and 
fully separable states are large enough. For all the factor considered, the only factor to increase the entanglement 
is the inter-mode amplification. 
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